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We discuss the phenomenology of initial-state parton-fc^ broadening in 
direct-photon production and related processes in hadron collisions. After 
a brief summary of the theoretical basis for a Gaussian-smearing approach, 
we present a systematic study of recent results on fixed-target and collider 
direct-photon production, using complementary data on diphoton and pion 
production to provide empirical guidance on the required amount of kx broad¬ 
ening. This approach provides a consistent description of the observed pattern 
of deviation of next-to-leading order QCD calculations relative to the direct- 
photon data, and accounts for the shape and normalization difference between 
fixed-order perturbative calculations and the data. We also discuss the uncer¬ 
tainties in this phenomenological approach, the implications of these results 
on the extraction of the gluon distribution of the nucleon, and the comparison 
of our findings to recent related work. 


1 


INTRODUCTION 


Direct-photon production has long been viewed as an ideal vehicle for measuring the 
gluon distribution in the proton [|l]. The quark-gluon Compton scattering subprocess 
(gq-’-'jq) dominates 7 production in all kinematic regions of pp scattering, as well as for 
low to moderate values of parton-momentum fraction, x, in pp scattering; the cross sections 
have been calculated to next-to-leading order (NLO) |Sj. The gluon distribution in the 
proton is relatively well-constrained at small x (x < 0 . 1 ) by deep-inelastic scattering (DIS) 
and Drcll-Yan (DY) data, but less so at larger x [[]]. Consequently, direct-photon data from 
fixed-target experiments that have been incorporated in several modern global parton distri¬ 
bution function analyses can, in principle, provide a major constraint on the gluon content 
at moderate to large x |||-(7j. 

However, a pattern of deviations between the measured direct-photon cross sections and 
NLO calculations has been observed 0. The discrepancy is particularly striking in the 
recently published higher-statistics data from E706 experiment E706 observed large 
deviations between NLO calculations and data, for both direct-photon and 7 r° inclusive 
cross sections. The final direct-photon results from UA 6 |T(J also exhibit evidence of similar, 
although smaller, discrepancies. The suspected origin of the disagreements is from effects of 
initial-state soft-gluon radiation. Such radiation generates transverse components of initial- 
state parton momenta, referred to in this discussion as kx- To be precise, as in [f|, kx 
denotes the magnitude of the effective transverse momentum vector, kx, of each of the two 
colliding partons. 

Evidence of significant kx has long been observed in measurements of dimuon, diphoton, 
and dijet pairs. A collection of measurements of the average transverse momentum of the 
pairs (( px) P air ) is presented in Fig. [1], for a wide range of center-of-mass energies (a/s) P~TH . 

The values of ( Px)pair are large, and they increase approximately logarithmically with 
increasing yf s . The values of ( k T ) per parton (estimated as ~(px) P air/ V 2) indicated by 
these DY, diphoton, and dijet data, as well as the inclusive direct-photon and 7 r° production 
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FIG. 1. (pt) of pairs of muons, photons, and jets produced in hadronic collisions versus yVs. 

data, are too large to be interpreted as “intrinsic” — i.e., due only to the finite size of the 
proton. (From these data, one can infer that the average kx per parton is about 1 GeV/c at 
fixed-target energies, increasing to 3-4 GeV/c at the Tevatron collider. One would expect 
( k T } values on the order of 0.3-0.5 GeV/c based solely on proton size.) Perturbative QCD 
(PQCD) corrections at NLO level are also insufficient to explain the size of the observed 
effects, and full resummation calculations are required to describe DY and W/Z jT2Hl4f, and 
diphoton [|15| , [I6| distributions. Values of (px)pair for DY and diphoton data exhibit similar 
trends versus energy; for DY data, pion-beam values are somewhat larger than those for 
proton beams at the same yVs. The dijet data hint at somewhat larger values of ( Px) P air 
than DY results at the same energy, a difference that may be related to different color-flow 
between initial and final states in DY and in dijet events, as well as to a larger contribution 
from gluon-induced subprocesses for dijet production. Similar soft-gluon (or k T ) effects are 
expected to be present in all hard-scattering processes, such as the inclusive production of 
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This paper presents a phenomenological model for kx effects in direct-photon production 
and, by extension, in all hard-scattering processes. Quantitative comparisons of this model 
with data from E706 have been reported previously 00- We will discuss the successes 
and uncertainties of this prescription, as well as the implications for determining the gluon 
distribution. 


hr SMEARING FORMALISM: THEORY AND PRACTICE 


We now briefly describe the theoretical underpinnings of kx effects using the Collins- 


Soper-Sterman (CSS) formalism |T9[. In this formalism, the px spectrum in hard-scattering 
processes is written as the convolution of the parton distributions /*//,., the C functions 
(representing the finite pieces of the virtual corrections), and two Sudakov form factors, 
S p and S NP . S p can be regarded as being perturbative in character and S NP as non- 
perturbative. The perturbative Sudakov form factor represents a formal resummation of 
soft-gluon emissions. The non-perturbative Sudakov form factor is determined from a fit 
to the data, is expected to be universal (for a given parton flavor and hadron type), and is 
usually parameterized as a Gaussian distribution. The dividing line between the perturbative 
and non-perturbative contributions is somewhat arbitrary (similar to the better known cases 
of parton distributions and fragmentation functions), and it is quantified by a theoretical 
scale in the resummation formalism. The px distribution (e.g., of the Droll-Yan pair) can 
be written symbolically as 

da 


dp T 


= cr o e 


~(S NP +S P ) 


{paii ® fi/hi'j (C a2 j <8> fj/hz) + ^ • 


( 1 ) 


The Y function is added to ensure a smooth matching between the low and the high-p^ 
regions, where the resummed and the fixed-order descriptions work better, respectively. 

At collider energies, most of the kx can be attributed to perturbative soft-gluon emission. 
However, for fixed-target kinematics, almost all of the kx is due to the non-perturbative 
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mechanisms. A proper treatment requires both the appropriate data to determine the non- 
perturbative input, and an implementation of the soft-gluon resummation formalism for the 
particular process. 


The resummation calculation for multiple soft-gluon emission in direct-photon produc¬ 
tion is quite challenging. The production rate and kinematic distributions of photon pairs 


produced in hadron interactions have already been calculated |16(], and a similar calculation 
of the transverse momentum distribution of a photon-jet system is also plausible, but more 
involved, since the final-state parton takes part in soft-gluon emission and in color exchange 
with the initial-state partons. A recent work on this subject |20| addressed only the effects 
of multiple soft-gluon radiation in the initial state. Incorporating jet definition in the for¬ 
malism is also not a fully resolved issue. Finally, the calculation of individual transverse 
momenta of the photon and the jet is further complicated by the fact that several over¬ 
lapping power-suppressed corrections can contribute. In the absence of the full resummed 
calculation, approximations are made in order to compare theory with data. 

In lieu of a rigorous calculation of the resummed px distribution, effects of soft-gluon 
radiation can be approximated by a convolution of the leading-order cross section with a 
/rr-smearing function ]2T[] . In the formalism described above, this is equivalent to absorbing 
all of the perturbative gluon emissions into the non-perturbative Sudakov form factor. Since 
no explicit resummation of soft-gluon emissions is performed, the average value of kx used in 
the smearing should be representative of the value observed (or expected) in the kinematic 
regime of the experiment. 

The expression for the leading-order (LO) cross section for direct-photon production at 
large px has the form: 
d 3 a 


E~ 


dp [ 


■{h 1 h 2 ^ 1 X) = 


ai a2(i3 


J2 [ dx l dx 2 f ai /hA x uQ 2 ) fa 2 /h 2 ( x 2,Q 2 )-^ (cqa 2 ^7%) S(s + t + u), 
. J vr dt 


( 2 ) 


where da/dt is the hard-scattering matrix element, and f ai /h-i an d f a2 /h 2 are the parton 
distribution functions (pdf) for the colliding partons a\ and a 2 in hadrons h\ and h 2 , re- 
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spectively. To introduce kx degrees of freedom, one extends each integral over the parton 
distribution functions to the kx- space, 


dx i f ai /hAxi,Q 2 ) -> dx 1 d 2 k Tl g(k Tl ) f ai /hAxi,Q 2 ), 


(3) 


(a corresponding substitution is done for parton a 2 in hadron h 2 ). 

The distribution g[kx) is usually taken to be a Gaussian, 

e -fer/( fc T> 

g(kT) = (4) 

where (k^) is the square of the 2-dimensional (2D) RMS width of the kx distribution for one 
parton (crf parton 2D ), and is related to the square of the 2D average of the absolute value of 
kx of one parton through (kj) = 4(/c r ) 2 /7r. We emphasize that ( kx) represents the average 
effective 2D transverse momentum per parton. (The average transverse momentum of the 
parton pair is, of course, a factor of y/2 larger than the average transverse momentum per 
parton.) 

The 4-vectors of the colliding partons are expressed in terms of the momentum fraction 
x of the partons. Ignoring parton and hadron masses, 

Xl = (E 1 +p h )/y/s, (5) 

where the parton four-vector is 

Pi = (E^kxnPh), (6) 


with 


and 




k T ! 

Xia/sJ ’ 


Ph 


1 

2 


Xi>/s 


k T i ' 

Xi y/s_ 


(Similar expressions are used for parton a 2 .) 


(7) 


( 8 ) 
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It is straightforward to evaluate the invariant cross sections, including kr effects, accord¬ 
ing to the above prescription. In general, because the unmodified PQCD cross sections fall 
rapidly with increasing px, the net effect of the kx smearing is to increase the yield. We 
denote the enhancement factor as K(pr)- Since the invariant cross section for direct-photon 
production is now a six-dimensional integral, it is convenient to employ Monte Carlo tech¬ 
niques in the evaluation of K(pr) ■ An exact treatment of the kinematics can be implemented 
in a Monte Carlo framework, but it is more difficult in an analytic approach. 

A Monte Carlo program that includes such a treatment of kx smearing, and the leading- 


order cross section for high-pr particle production, has long been available |2l[. The program 
provides calculations of many experimental observables, in addition to inclusive cross sec¬ 
tions. The program can be used for direct photons, jets, and for single high-p T particles 
resulting from jet fragmentation (such as inclusive 7 r° production). Unfortunately, no such 
program is available for NLO calculations, but one can approximate the effect of kx smearing 
by multiplying the NLO cross sections by the LO /^-enhancement factor. Admittedly, this 
procedure involves a risk of double-counting since some of the /^-enhancement may already 
be contained in the NLO calculation. However, such double-counting effects are expected to 
be small, and consequently this uncertainty on K(px ) should also be small. (For example, 
the NLO estimate for ( px) P air of the diphoton pairs produced in pp —> 77 at y/s = 31.5 
GeV is on the order of a few hundred MeV/c, while the resummed prediction is well over 1 
GeV/c 0.) 

It is clear that this type of treatment of kx effects is model dependent. In particular, 
different functional forms can lead to quantitatively different answers. For example, adding 
substantial non-Gaussian tails in kx smearing can affect the output distributions. One of 
the strengths of the approach we follow in this paper is that the kx distribution used in the 
smearing is based on experimental information, and the Gaussian character of the kx effects 
is consistent with the data observed by E706 [Q. Moreover, any non-Gaussian tails may 
result primarily from single hard-gluon emission, and such contributions should therefore 
already be contained in the NLO cross sections used in this analysis. 
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A complete treatment of soft-gluon radiation in high-pr production, including the ap¬ 
propriate non-perturbative input, should eventually predict the effective k T values expected 
for each process and y/s. We will employ (hr) values representative of the kinematic distri¬ 
butions in the data, and based on comparisons with the same model as used to modify the 
NLO inclusive cross sections. 

The effects of soft-gluon radiation are also included in QCD Monte Carlo programs such 
as Pythia [^2|, Isajet f23|, and Herwig |24|]. However, in these programs the emission is 
normally cut off at a relatively high parton virtuality, with the remaining kx effect supplied 
by a Gaussian smearing similar to that discussed above. For fixed-target energies, essentially 
all of the kx effects are provided by this phenomenological Gaussian term. The above 
programs differ in the details of the way parton energy and momentum are rescaled after 
k T is inserted, which can also produce quantitative differences in results. 


APPLICATIONS OF THE k T MODEL TO DATA 

The experimental consequences of kx smearing are expected to depend on the collision 
energy. At the Tevatron collider, the smallest photon px values probed by the CDF and 
D0 experiments are rather large (10-15 GeV/c), and the ^-enhancement factors modify 
only the very lowest end of the px spectrum, where px is not significantly larger than 
kx • In the energy range of the E706 measurements, large /^-effects can modify both the 
normalizations and the shapes of the cross sections as functions of px■ Consequently, E706 
data provide a particularly sensitive test of the kx model. At lower fixed-target energies, 
the kx enhancements are expected to have less px dependence over the range of available 
measurements, and can therefore be masked more easily by uncertainties in experimental 
normalizations and/or choices of theoretical scales. Nonetheless, the UA6 and WA70 data 
generally support the expectations from kx smearing. 






Comparisons to Tevatron Collider Data 


At the Tevatron collider, the above model of soft-gluon radiation leads to a relatively 
simple modification of the NLO cross section. In Fig. we compare the CDF and D0 
isolated direct-photon cross sections [^5|] to theoretical NLO calculations with and without kx 
enhancement. In the lower part of the plot we display the quantity (Data-Theory)/Theory; 
for the collider regime we did not multiply the NLO theory by the /^-enhancement factor, 
but instead displayed the full deviation of the NLO calculation from the data. The expected 
effect from kx enhancement is also shown for (kx) = 3.5 GeV/c. This is the approximate 
value of ( kx ) per parton measured in diphoton production at the Tevatron [|5|], and one 
expects a similar (kx) per parton for single-photon production. (In the diphoton process, 
the 4-vectors of the photons can be measured precisely, providing a direct determination of 
the transverse momentum of the diphoton system, and thereby (kx)-) 

As seen in Fig. ^ the kx effect diminishes rapidly with px and is essentially negligible 
above ~30 GeV/c. The trend of deviations of NLO calculations from the measured inclusive 
cross sections is described reasonably well by the expected kx effect. Some of the observed 
excess can be attributed to the fragmentation contribution to isolated direct-photon pro¬ 
duction |26j] , but this alone cannot account for the entire deviation of the theory from data. 

The larger statistics in the Tevatron collider Run IB samples (currently under analysis) 
will allow for a more detailed examination of the \ow-px behavior of the photon cross sec¬ 
tion. (The CDF data included in the plot are from Run IA only.) A similar enhancement 
is expected for jet production at low px, but larger experimental uncertainties, and the 
relatively large additional non-perturbative effects expected in this region [|27|, preclude a 
useful comparison. 
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FIG. 2. Top: The CDF and D0 isolated direct-photon cross sections, compared to NLO theory 
without hr (dashed) and with hr enhancement for (hr) = 3.5 GeV/c (solid), as a function of pr- 
Bottom: The quantity (Data-Theory)/Theory (for theory without kr adjustment), overlaid with 
the expected effect from kr enhancement for (kr) = 3.5 GeV/c. The error bars have experimental 
statistical and systematic uncertainties added in quadrature. 
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Comparisons to E706 Data 


The conventional NLO calculations yield cross sections that are signficantly below the 
E706 direct-photon and tt° measurements (see Figs. |. f|, andj^). No choices of current 
parton distributions, or conventional PQCD scales provide an adequate description of the 
data (for the comparisons presented here all QCD scales have been set to px/ 2). The 
previously described /cr-enhancement algorithm was used to incorporate the effects of soft- 
gluon radiation in the calculated yields. That is, the theory results plotted in the figures 
represent the NLO calculations multiplied by /^-enhancement factors K{px ) PB| . 

Because parton distributions for nucleons are known better than those for pions, we first 
present comparisons of the various model calculations with proton-beam data. As seen at 
the bottoms of Figs. [| and [I[ the NLO theory, when supplemented with appropriate k T 
enhancements, is successful in describing both the shape and normalization of the E706 
direct-photon cross sections at both y/s = 31.6 GeV and 38.8 GeV. As expected, the kx- 
enhancement factors affect the normalization of the cross sections, as well as the shapes of 
the px distributions. The values of ( kx) = 1.2 GeV/c at yfs = 31.6 GeV, and 1.3 GeV/c at 
a/s = 38.8 GeV, provide good representations of the incident-proton data. Both kx values 
are consistent with those emerging from a comparison of the same PQCD Monte Carlo with 
E706 data on the production of high-mass 7 r° 7 r°, 77 r° and 77 pairs 0. 

In Fig. |H we present comparisons between calculations and E706 data for 7 r”Be interac¬ 
tions at yfs = 31.1 GeV. Here again, the uncorrected {{kx) = 0) NLO theory is not consistent 
with the data. Once the /c^-enhancement factors are applied, good agreement is observed 
between data and calculations for {kx)~ 1.4 GeV/c. (Note that DY data lead one to expect 
a higher value of {kx) for a n~ beam than for a proton beam of the same energy.) 

For comparison, results of calculations using {kx) values ±0.2 GeV/c of the central values 
are also shown in the figures. These can be taken as an indication of uncertainties on {kx) 
(see next section). The corresponding enhancement factors K{px ) at yfs = 31.6 GeV are 
displayed in Fig. || 
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FIG. 3. Top: The photon and 7r° cross sections from E706 at \fs = 31.6 GeV compared to 
Ayr-enhanced NLO calculations. Bottom: The quantity (Data- Theory)/Theory for direct-photon 
production, using Ayr-enhanced NLO calculations for several values of (Ayr)- The error bars have ex¬ 
perimental statistical and systematic uncertainties added in quadrature. The points corresponding 
to calculations with different (kx) are slightly staggered in px , to reduce the overlap of experimental 
error bars. 
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FIG. 4. Top: The photon and 7r° cross sections from E706 at \fs = 38.8 GeV compared to 
/cT’-enhanced NLO calculations. Bottom: The quantity (Data- Theory)/Theory for direct-photon 
production, using /^-enhanced NLO calculations for several values of {hr) ■ The error bars have 
experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 5. Top: The photon and 7r° cross sections from E706 at \fs = 31.1 GeV for incident 7T~ 
beam, compared to /^-enhanced NLO calculations. Bottom: The quantity (Data-Theory)/Theory 
for direct-photon production, using fey-enhanced NLO calculations for several values of {kx)- The 
error bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 6. The variation of kx enhancements, I\(px), relevant to the E706 direct-photon data for 
protons at = 31.6 GeV, for different values of (kx)- 

It is interesting to note that, in this energy range, the kx enhancement does not decrease 
with increasing px (as for the case of the collider), and, in fact, increases at the highest 
values of px- This is a consequence of the wide x-range spanned by E706, and can be 
understood through the following argument. At lowpr, a (kx) of 1.2 GeV/c is non-negligible 
in comparison to the px in the hard-scattering, and the addition of kx smearing therefore 
increases the size of the cross section (and steepens the slope). At high px (corresponding 
to large x ), the unmodified NLO cross section becomes increasingly steep (due to the rapid 
fall in parton densities), and hence the effect of kx smearing again becomes larger. 

NLO calculations for 7r° production have a greater theoretical uncertainty than those for 
direct-photon cross sections since they involve parton fragmentation. However, the kx effects 
in 7r° production can be expected to be generally similar to those observed in direct-photon 
production, and the 7r° data can be used to extend tests of the consequences of kx smearing. 
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Figures |3|, f|, and 0 also show comparisons between NLO calculations |[29]| and 7r° production 
from E706, using BKK fragmentation functions (ff) |30|. The previously described Monte 
Carlo program was employed to generate A/r-enhancement factors for 7r° cross sections, and 
(kx) per parton values similar to those that resulted in good agreement with direct-photon 
data also provide a reasonable description of 7T° production. For 7T° production, an additional 
smearing of the transverse momentum, expected from jet fragmentation, has been taken into 
account. 


Comparisons to WA70 and UA6 Data 


We have examined the expectations for the size of soft-gluon effects for fixed-target ex¬ 
periments WA70 and UA6. Both experiments have measured direct-photon production with 
good statistics, and their data have been included in recent global fits to parton distribu¬ 
tions. WA70 measured direct-photon and 7r° production in pp and n~p collisions at a/s 
= 23.0 GeV |31|], and UA6 has recently released [HJ their final results (with substantially 
reduced uncertainties) for direct-photon production in pp and pp collisions at a/s = 24.3 
GeV. These center of mass energies are somewhat smaller than those of E706, and the (kx) 
values are therefore expected to be smaller (perhaps of the order of 0.7-0.9 GeV/c, based 
on Fig. 0). WA70 has compared kinematic distributions observed in diphoton events (for 
7 r~p interactions) to NLO predictions, and has found that smearing the NLO theory with 
an additional ( kx) of 0.9±0.2 GeV/c provides agreement with their data [^]. We therefore 
use this (kx) as the central value for the /cr-enhancement factors for both experiments, and 
vary the (kx) by ±0.2 GeV/c, as with E706. The corresponding kx enhancement expected 
over the range of measurements is shown in Fig. |7j. Over this narrower px range, the effect 
of kx is essentially to produce a shift in normalization. 

Comparisons of the WA70 direct-photon and 7r° cross sections with the fc^-enhanced 
NLO calculations are shown in Figs. ^ and 0. Renormalization and factorization scales of 
Px/ 2 are used in the NLO calculations, as in the E706 comparisons. The 7r° cross sections 
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FIG. 7. The A^-enhancement factors for direct-photon production expected for WA70 pp data. 

both for incident proton and beams, and the photon data from incident beam, all lie 
above the NLO calculations for (kx) — 0, and are in better agreement with the A^-enhanced 
calculations; only the photon cross section for incident protons seems not to require a kx 
correction. 

The latest photon cross sections from UA6 for pp and pp scattering are shown in Figs. E3 
and [id]. The photon cross section for pp interactions lies clearly above the NLO calculation 
for (k T ) = 0, but is consistent with /c T -adjusted calculations for (k T ) in the range of 0.7-0.9 
GeV/c. The result for pp interactions is also above the unmodified NLO calculation, but 
requires a smaller value of (kx)- We note that the dominant production mechanisms for 
the two processes are different: quark-gluon Compton scattering dominates for pp, and qq 
annihilation for pp at the UA6 energy. As in the case of E706 and WA70 measurements, 
the UA6 7r° cross sections are also higher than the NLO calculation without kx-, and can be 
much better described by introducing kx enhancement. 
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FIG. 8. Top: The photon and 7r° cross sections from WA70 at y / i = 23.0 GeV for incident pro¬ 
tons compared to /^-enhanced NLO calculations. Bottom: The quantity (Data-Theory)/Theory 
for direct-photon production, using fey-enhanced NLO calculations for several values of (kx)- The 
error bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 9. Top: The photon and 7r° cross sections from WA70 at \/s = 23.0 GeV for incident 7r~ 
beam, compared to /^-enhanced NLO calculations. Bottom: The quantity (Data-Theory)/Theory 
for direct-photon production, using fey-enhanced NLO calculations for several values of {kx)- The 
error bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 10. Top: The photon and 7r° cross sections from UA6 at y''s = 24.3 GeV for an 
incident proton beam compared to for-enhanced NLO calculations. Bottom: The quantity 
(Data- Theory)/Theory for direct-photon production, using fop-enhanced NLO calculations for sev¬ 
eral values of (fop)- The error bars have experimental statistical and systematic uncertainties added 
in quadrature. 
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cident antiproton beam compared to fc^-enhanced NLO calculations. Bottom: The quantity 
(Data- Theory)/Theory for direct-photon production, using Ly-enhanced NLO calculations for sev¬ 
eral values of (fcy). The error bars have experimental statistical and systematic uncertainties added 
in quadrature. 
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DISCUSSION OF /ct-SMEARING PROCEDURES 


In this section, we consider uncertainties in the Monte Carlo method employed in our 
previous discussion, and comment on an analytic approach for calculating such kx enhance¬ 
ments. 


Uncertainties in the Monte Carlo Model and Related Issues 

Our approach contains several phenomenological parameters that affect the range of 
its results, the most important of which is the amount of Gaussian smearing represented 
by (kx)- Thus far, both the value of (k T ) (with its dependence on y/s) and its possible 
range of uncertainty can only be determined empirically. A range of variation of ±0.2 
GeV/c for (kx)~l—2 GeV/c appears reasonable on the basis of several considerations. These 
include: (i) the range of kx values inferred from the E706 high-mass pair distributions for 
several kinematic variables; (ii) observed differences between photon and 7r° results; (iii) 
comparisons of (kx) values required in inclusive cross sections to those representing the 
properties of massive pairs at E706 and WA70/UA6 energies, and (iv) differences between 
dimuon, diphoton, and dijet values of ( Px) P air • (The dependence of the /^-enhancement 
factor on the variation in (kx) was illustrated in Figs. ^ and [7].) 

For the fixed-target experiments discussed in the previous section, we already presented 
calculations using (kx) values of 0.2 GeV/c above and below the selected central values. 
The observed variation in predictions reflects an uncertainty in the fey-enhanced theoretical 
results. The dependence of K(px) on (kx) is especially strong at the extremes (low and high) 
in px- These calculations used renormalization and factorization scales of px/2- Changing 
the scale to px reduces predicted cross sections at fixed-target energies by about 30-40%, a 
factor comparable to the estimated spread in the kx factors over much of the px range; the 
full uncertainty in the calculations must include contributions from uncertainties in both 
(kx) and QCD scales. 
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Another parameter in the model, a propagator mass m, is introduced to regularize the 
divergences of the leading-order QCD matrix elements due to the propagator factors 1/s, 
1/t, and 1 ju in the configuration where one or more of these invariants approaches zero. To 
avoid the large weights associated with these configurations, the propagators are replaced 
by l/(s + m 2 ), 1 /(t - m 2 ), and 1 / (u - m 2 ), respectively, where m has a typical value of order 
1 GeV. The physical effect of the propagator mass is to cut off the region where almost all 
of the transverse momentum of the produced photon is due to the Gaussian smearing, and 
very little to the hard scattering. While the /^-enhancement factor is relatively insensitive 
to the value of m for the measured range at collider energies, it is somewhat sensitive for 
data at fixed-target energies. 

To illustrate the sensitivity of the /^-enhancement factor K(px) to the choice of the 
propagator mass, we display in Fig. |1^ the results for m = 1.0 GeV (default value), 1.3 
GeV, and 0, for pBe—at a/s = 31.6 GeV. Above px values of 5.5-6.0 GeV/c, the kx 
enhancement has little sensitivity to m. At lower px, the dependence on the propagator mass 
should be considered as a measure of an uncertainty of the model in this region of phase 
space. Clearly, larger enhancements are obtained at low px when there is no propagator 
mass (m — 0). 

The value of (kx), appropriate in the calculation, depends on the kinematic situation. 
The increase of the value of (kx) with s is understood, in general terms, as the result of 
an increase in the available phase space for multiple-soft-gluon emission. This growth is 
predicted by a CSS-type of resummation for Drell-Yan processes |33j and, as illustrated 
in Fig. |I|, has been observed in Drell-Yan and diphoton data. For simplicity, the model 
calculation assumes a constant value of (kx) for a given y/s. In reality, various physical effects 
can modify this simple choice and cause a modification in the shape of the enhancement 
factor K(px). Below, we discuss a few of these effects. 

First, one might expect a dependence of (kx) on s, the hardness of the partonic inter¬ 
action, similar to that on s. A study of diphoton production at y/s = 31.5 GeV, using the 


RESBOS program |Tj], indicates that the (kx) increases from about 1.2 GeV/c for photons 
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FIG. 12. The effect of different propagator masses m for predictions of kx enhancements K(px) 
for E706 direct-photon data at yfs = 31.6 GeV. The kx enhancement used by the MRST group 
fill is also shown in the figure. 

with px of 5 GeV/c to approximately 1.4 GeV/c for photons with px of 10 GeV jc. If direct- 
photon production were to have a similar behavior, the anticipated kx enhancement at a px 
of 10 GeV/c would be about 40% higher than for a constant (kx) value of 1.2 GeV/c (see 
the dependence of the enhancement factor on (kx) in Fig. |5|). 

Suppression mechanisms may also exist at large px due to the restriction of phase space 
for gluon emission from large-x partons. A gluon emitted from an incoming parton carries 
away a momentum fraction Ax = 2pxe y /y/s, where y is the rapidity of the gluon. Since 
gluons emitted at forward rapidities would carry away more of the incident parton’s mo¬ 
mentum than available, the effective rapidity range for gluon emission becomes restricted to 
central rapidities. Such effects are not taken into account correctly in the simple kx models 
discussed above. 
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Finally, there has been much recent interest in studying the effects of resumming large 
logarithms of the form ln(l —Xt), where Xt = 2p T /y/s |34|j35|l . As Xx approaches 1, for any 
hard-scattering process, the perturbative cross section is enhanced by powers of ln(l — xt) 
that have to be resummed at all orders. These types of effects should currently be negligible 
for direct-photon production at the Tevatron collider (because data are available only for 
relatively small values of xt), but may be important at fixed-target energies. Definitive 
answers to the question of whether this effect leads to a net enhancement or suppression, 
and its consequences for the shape of the px distribution, are not presently available. A 
treatment that includes both kx and threshold resummation effects may be necessary for a 
more satisfactory description of the fixed-target data. 


Analytic Smearing Methods 


An alternative to the Monte Carlo calculation of the kx enhancement is its implemen¬ 
tation through analytic convolution of the theoretical cross section with a (Gaussian) kx 
distribution in either one or two dimensions. The latter is a better approximation, but 
the correction to the one-dimensional convolution is second order in kx/px , leading to a 
difference of about 10% for the fixed-target applications. To compare to the Monte Carlo 
results, the parameters of such analytic convolutions must be interpreted in terms of the 
parton-level (kx) values. 

We consider the kinematics of the production of a direct photon accompanied by a recoil 
jet. As mentioned before, the total transverse momentum imparted to the final state ( 7 +jet) 
by the colliding partons has an RMS width 


&2partons ,2D 


V2ui p 


arton,2D • 


(9) 


However, the single photon is subject only to half of this transverse kick, because it shares 
the total px with the recoil jet. Consequently, the kx RMS width relevant for an analytic 
smearing of the single-inclusive photon cross section is 


25 
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(10) 


The above can be illustrated by a simple example. Consider the production of a photon 
and a jet with equal and opposite px values (for example, 4 GeV/c) close to 90 degrees in 
the center of mass. Compare this to the situation when the colliding partons impart some 
amount of kx, say 1 GeV/c, in the direction of the photon. The photon px is now 4.5 
GeV/c, and the jet px becomes 3.5 GeV/c (in the opposite direction), resulting from a total 
of 1 GeV/c of hr imparted to the photon+jet system. Thus, the photon itself only receives 
half of the partonic total. This factor of two is critical, and can be easily overlooked, in 
calculating kx smearing using analytic methods. The above conclusions have been verified 
within the explicit Monte Carlo treatment of event kinematics. (The convolution formulae 
for kx smearing are discussed more fully in the Appendix.) 

Recently, the MRST group has produced a new set of parton distributions, incorporating 
kx effects in their analysis of the WA70 and E706 direct-photon data |I6|; they used an 
analytic smearing technique rather than a Monte Carlo approach. The correction algorithm 
differs in detail with the one discussed above, and seems to require a significantly smaller 
( kx ) to describe the data (about 0.65 GeV/c per incoming parton as compared to 1.2 GeV/c 
that was used for E706 at a/s = 31.6 GeV). This discrepancy may be due to a difference 
in the interpretation of kx that is related to the factor of two present in Eq. [H]. If we 
reinterpret the ( kx )mrst of 0.65 GeV/c as our 1.3 GeV/c per parton, then agreement is 
restored. Nevertheless, the /or-enhancement factors differ somewhat in the two techniques, 
leading to different conclusions about the gluon distribution, especially at large x (see next 
section). For comparison, the MRST kx factor for E706 (at a/s = 31.6 GeV) has also been 
plotted in Fig. [L2]; it is similar in size to our calculations for ( kx ) values in the range of 
1.2-1.4 GeV/c, but has a different shape in px■ Until these differences are understood, we 
must assume that there is a significant shape-dependent uncertainty associated with the 
particular assumptions used for modeling kx effects. 
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IMPACT ON THE GLUON DISTRIBUTION 


It is now generally accepted that the uncertainty in the gluon distribution at large x 
is still quite large. Thus, it would appear important to incorporate further constraints 
on the gluon, especially from direct-photon data. In the following, we describe a global 
pdf fit that employs our kx model in the analysis of E706 direct-photon measurements. 
Since jet production at the Tevatron collider is another available constraint for the gluon 
content at large x, we discuss the consistency between the E706 direct-photon results and 
the CTEQ4HJ gluon distribution (derived using high-p T jet data from CDF) p7|. 


Application of kx Enhancements in a pdf Fit 


To investigate the impact of kx effects on determinations of the gluon distribution, we 
have included the E706 direct-photon cross sections for incident protons, along with the 
DIS and DY data that were used in determining the CTEQ4M pdfs, in a global fit to the 
parton distribution functions. The CTEQ fitting program was employed to obtain these 


results |38|], using the NLO PQCD calculations for direct-photon cross sections, adjusted by 
the /^-enhancement factors. However, the WA70, UA6, CDF, and D0 data were excluded 
from this particular fit. The resulting gluon distribution, shown in Fig. [H| is similar to 
CTEQ4M, as might have been expected, since the /^-enhanced NLO cross sections using 
CTEQ4M provide a reasonable description of the data shown in Figs. and |4]. 

The data sets used in determining CTEQ4M did not include the E706 direct-photon 
cross sections, but did use earlier direct-photon data from UA6 and from WA70 (without 
kx corrections), along with the inclusive jet cross sections from CDF and D0 |h|. The jet 
cross sections were particularly useful for defining the gluon distribution in CTEQ4M at 
moderate values of x. 

The new MRST gluon distribution (also shown in Fig. [13]) is significantly lower than 
CTEQ4M (and MRSR2) at large x. While the MRST fit employs kx enhancements, it 
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FIG. 13. A comparison of the CTEQ4M, MRST, and CTEQ4HJ gluons, and the gluon dis¬ 
tribution derived from fits that use E706 data. The g] and g{ gluon densities correspond to the 
maximum variation in (kx) that MRST allowed in their fits. 

attempts to accommodate the WA70 incident-proton direct-photon data, which does not 
exhibit an obvious kx effect. In addition, the MRST /^-enhancement factors are larger 
at large px than our corresponding results, resulting in a smaller gluon at large x. This 
further serves to illustrate the extent to which the extracted gluon distribution is affected 
by the specific procedures applied in the fit. In contrast, the CTEQ4HJ specialized gluon 
distribution (discussed in more detail in the next subsection) is much larger than CTEQ4M 
in the same x range. The current spread of the solutions for the gluon distribution at large 
x is uncomfortably large, and additional effort is required to resolve these discrepancies. 







Discussion of the CTEQ4HJ Gluon 


As presented above, when analyzed with /^-enhancement factors, the E706 direct-photon 
data lead to a gluon distribution similar to that in CTEQ4M. The implications for the size 
of the gluon at large x are especially important because of the excess observed by CDF in 
the high-pT’ inclusive jet cross section (when compared to calculations using conventional 
pdfs). The CTEQ collaboration produced a global fit (CTEQ4HJ) j37| to improve the 
description of the high-pr jet data from CDF in Run IA [flTJ. The high-pr data points were 
given an enhanced weight to emphasize them in that fit. The resulting CTEQ4HJ parton 
distributions produce a jet cross section that by design follows the CDF data points more 
closely than the cross section obtained using CTEQ4M. 

The CDF inclusive jet cross section from Run IB |4^] demonstrates an excess at high p 7 - 
similar to that observed in Run IA. For D0, the inclusive jet cross section from Run IA+IB 
is consistent with the NLO QCD calculation using conventional parton distributions such 
as CTEQ4M, but can also be well-described with calculations using the CTEQ4HJ parton 
distribution functions; in fact, calculations using CTEQ4HJ result in better y 2 agreement 
with the D0 jet data J43] . 


The CTEQ4HJ quark distributions are similar to those obtained in a more standard fit 
(as for example CTEQ4M), since the DIS and DY data provide significant constraints on 
the quark distributions over all x (and on the gluon distribution at small x) ||. However, 
the CTEQ4HJ gluon distribution is significantly larger at high x than that of CTEQ4M, by 
a factor of 1.5 at x — 0.5 (for Q — 5 GeV/c), and by a factor of 5 at x = 0.7. 

Because of the dominance of the qq scattering subprocess in the Tevatron jet cross sections 
at high px, a large change in the gluon distribution is required to generate a relatively small 
change in the jet cross section. As the CTEQ exercise has demonstrated, until the theoretical 
issues related to interpretation of direct-photon data are resolved, there is freedom within the 
data sets used in the global fits to change the gluon distribution in this way. It should also be 
noted that a recent analysis |IT| of deuteron and proton structure functions, using corrections 
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FIG. 14. Top: The photon cross section from E706 at y/s = 31.6 GeV compared to ^-enhanced 
NLO calculations using the CTEQ4HJ parton distribution functions. Bottom: The quantity 
(Data- Theory)/Theory using /^-enhanced NLO calculations for several values of (for)- The error 
bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 15. Top: The photon cross section from E706 at yfs = 38.8 GeV compared to /^-enhanced 
NLO calculations using the CTEQ4HJ parton distribution functions. Bottom: The quantity 
(Data- Theory)/Theory using /^-enhanced NLO calculations for several values of (for)- The error 
bars have experimental statistical and systematic uncertainties added in quadrature. 


31 






for nuclear binding effects in the deuteron, suggests that the down-quark distribution in the 
nucleon at large x may be significantly larger than previously assumed; clearly this also 
influences the calculated cross sections for the high -px jet production. 


The CTEQ4HJ gluon distribution is compared to those from CTEQ4M, MRST, and the 
fits including the E706 data in Fig. [T3]. Figures [T4] and |T5] show comparisons of the E706 
direct-photon cross sections and the NLO calculations using the CTEQ4HJ parton distribu¬ 
tions, with and without the /^-enhancement factor. The shape in px of these calculations 
appears less consistent with the data than the corresponding results using the CTEQ4M 
gluon. However, current uncertainties in the understanding of hr effects in direct-photon 
production (discussed in the previous section), preclude an unambiguous interpretation of 
this difference. 


CONCLUSIONS 

We have described a phenomenological model for kx effects in which (kx) values used 
in the calculations of /^-enhancement factors are derived from data. Despite uncertainties, 
the results are remarkably successful in reconciling the data and theoretical calculations for 
a broad range of energies. The /^-enhancement factors improve the agreement of QCD 
calculations with E706, UA6, and n~ beam WA70 direct-photon cross sections over the full 
px range of measurements, and at the \ow-px end of CDF and D0 results. All fixed-target 
7 r° measurements also agree much better with such /^-enhanced calculations. 

The proper treatment of soft-gluon radiation in direct-photon cross sections can affect 
the extraction of the gluon distribution, especially at large values of x. In this particular 
treatment of kx enhancement, the E706 data, which span the widest x-range of any direct- 
photon experiment, are in better agreement with the gluon distribution from CTEQ4M than 
from CTEQ4HJ. However, within the phenomenological approach, any physical mechanism 
which gives rise to less enhancement at large x than our specific model calculation can make 
the CTEQ4HJ gluon more consistent with the E706 data. 
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A definitive conclusion regarding the quantitative role of hr effects in hard scattering, 
and reliable additional information on the large-x gluon distribution, awaits more complete 
theoretical calculations. The new generation of direct-photon and inclusive jet measure¬ 
ments serves as a strong impetus, and provides a testing ground, for new theoretical de¬ 
velopments, perhaps incorporating both kr and large -x resummation, which will have less 
model dependence than the formulations discussed in this paper. Further progress built on 
this interplay between theory and experiment will allow a more definitive determination of 
the gluon distribution, especially in the large-x region, where significant uncertainties still 
remain. 
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Appendix 

In this appendix, we collect formulae relevant to the analytic treatment of kj- effects in 
inclusive cross sections. As discussed in the main text, we assume a Gaussian description 
for the parton /^-distributions; the width parameters entering the formulae are labeled 
explicitly to help keep track of their meaning. 

For definiteness, let us consider direct-photon production. The full 2-dimensional con¬ 
volution of the (parametrized) differential cross section E (for example, E = da/dpx ) with 
the Gaussian fc-r-smearing functions can be written as: 

£' (p T ) = / d 2 k Tl d 2 k T2 d 2 q T e^V^) 1 e ~ k T 2 /( k T 2 ) 

J 7t(/4i) ^\ k T 2 ) 

xE (q T ) h (2) (Pr — qT ~ + ^t 2 )), (H) 

or, integrating out the (5-function constraints, as 

£'(pr) = [ d 2 k T — l —exp {-kl/a 2 2D ) E(|p r - k T |), (12) 

J 7T<J 'y,2D 

where ct^^d is the width parameter appropriate for the smearing of the direct-photon inclu¬ 
sive distribution; it is related to the width of the parton /^-distribution through Eq. (S3). 

The /^-enhancement factors, defined as K(p T ) = E'(p T )/E(p T ), can be calculated nu¬ 
merically. For a more intuitive treatment, one can simplify the discussion by reducing the 
above equation to a 1-dimensional case. Let us decompose the kx vector into components 
parallel and perpendicular to px- Then only the parallel component strongly affects the value 
of the photon px, while the perpendicular component affects px much less for typical kx^Pr 
configurations (since it adds to px in quadrature). Neglecting the effect of the perpendicular 
component, and integrating it out, one arrives at a 1-dimensional approximation: 
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where 
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(14) 


and we used Eq. ([TT]) to relate <t 7i2 d and (Ji par ton, 2 D- 

Two particularly interesting special cases are the exponential form of the cross section, 
used to describe the data at fixed-target energies, and the l/p? form, appropriate at colliders. 
For an exponential representation of the cross section, E ~ exp(— bpr), one obtains 


'■*- 1 , 2,2 \_ 1 u2, 


K(p T ) = exp(-6 P lri , r , mSD ) = ex P(^ b ( k r) )■ 


(15) 


Using the LO calculation of pp —> jX at y/s = 31.6 GeV (and {hr) = 0), one finds b ~ 
1.8 (GeV/c) -2 at pr = 6.5 GeV/c (a value in the middle of the p T range of E706). The 
approximate 1-dimensional formula above yields K & 2.1 for (hr) = 1.2 GeV/c, compared 
to K LO ~ 2.0 obtained in the full LO Monte Carlo calculation. In general, despite the 
approximations in the treatment, the analytical results are quite close to the results of the 
LO calculations. The factor in the relation of cx 7i 2 d and cri par ton, 2 D (Eq. fh]) is crucial to 
obtaining this consistency. 

A different representation, useful, for example, for parametrizing CDF and DO measure¬ 
ments, assumes E ~ 1 jpVp. For this parametrization (or more general functional forms) one 
can expand E(pt — kr) as a power series in kr (for kr small compared to pr)- 

E(p T - k T ) = S (p T ) + T^k^E\"{p T ) + -^E (4) (p T ) + ••• (16) 


(the odd powers of kr integrate out to zero). One obtains: 


K(p T ) 


1 _l ( k T ) 2 n{n + 1 ) 
27T pj, 


+ 


( k T ) 4 n(n + l)(n + 2 ){n + 3) 

87T 2 Py 


+ ... 


(17) 


For a constant (or a slowly changing) slope parameter n (and for {kr) Pr), the effects of 
k T smearing decrease as 1/p^, as might be expected for a power-suppressed process. 
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